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AN APPROXIMATE METHOD FOR SOLVING THE INVERSE 
SCATTERING PROBLEM WITH FIXED-ENERGY DATA 

'q^ • A. G. RAMM AND W. SCHEID 

Abstract. Assume that the potential q{r), r > 0, is known for r > a > 0, 
and the phase shifts Si{k) are known at a fixed energy, that is at a fixed k, 
^ ' for ( = 0, 1, 2, . . . . The inverse scattering problem is: find q{r) on the interval 

Q ' 0< r < a, given the above data. A very simple approximate numerical method 

, is proposed for solving this inverse problem. The method consists in reduction 

■ of this problem to a moment problem for q{r) on the interval r £ [0, a]. This 

00 ' moment problem can be solved numerically. 

(N 

> 

in 

I 1. Introduction 

Finding a potential q{r), r = \x\, from the phase shifts 6i{k) for angular quantum 
numbers I = 0, 1, 2, ... , known at a fixed wave number k > 0, that is, at a fixed 
0^ . energy, is of interest in many areas of physics and engineering. A parameter-fitting 

' procedure for solving this problem was proposed in the early sixties by R. G. Newton 

^ , and discussed in This procedure has principal drawbacks which have been 

pL|' discussed in detail in the paper [|l|. 

In [1] AGR pointed out for the first time that the Newton - Sabatier method 
widely used by the physicists to invert fixed-energy phase shifts for the potential 
is not really an inversion method but a parameter-fitting procedure; that this pro- 
cedure is not always applicable to the data; and that it is not proved that this 
procedure leads to the original potential, which produced the original phase shifts. 
Moreover, the Newton - Sabatier procedure, if it is applicable, can only produce a 
^ . , potential which is analytic in a neighborhood of (0, -|-cx)), with a possible pole at 

' r = 0. A general numerical method is given in [1] to construct piecewise-constant 

compactly supported spherically-symmetric potentials which are quite different but 
have nearly the same sets of phase shifts at a fixed energy. 

If A; > is large, then one can use a stable numerical inversion of the fixed-energy 
scattering data in the Born approximation ]To| , Section 5.4]. Error estimates for the 
Born inversion are obtained in [To[ |. 

An exact (mathematically rigorous) inversion method for solving a 3D inverse 
scattering problem with fixed-energy noisy data is developed in |ll[ where error 
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estimates were derived and the stability of the solution with respect to small per- 
turbations of the data was estimated. This method, although rigorous, is not very 
simple. 

The aim of this paper is to propose a novel, approximate, quite simple in prin- 
ciple, method for inverting the data {5i}, Z = 0, 1, 2, . . . , given at a fixed fc > 0, for 
the potential q{r). In most physical problems one may assume that q{r) is known 
for r > a (near infinity), and one wants to find q{r) on the interval [0,o], where 
a > is some known radius. 

Our basic idea is quite simple: since 6i and q(r) for r > a are known, one can 
easily compute the physical wave function tjji (r) for r > a, and so the data 

(1.1) {VKa),V'Ka)}, i = 0,1,2,... 

can be obtained in a stable and numerically efficient way (described in Section 2 
below) from the original data {Si}, {q{r) for r > a}. 

Now we want to find q{r) on the interval r £ [0, a] from the data (1.1). The 
function tpi{r) on the interval [0,oo) solves the integral equation 

(1.2) iPiir)=^l°\r)- f gi{r,p)q{p)Mp)dp, 0<r<a 

Jo 

where gi is defined in (A. 12) and tpl°\r) is a known function which is written 
explicitly in formula (2.7) below. Let 

(1.3) -iMo) - 4'\a)] bi, -[i4{a) - i^f^'{a)] -.^ pu ^ = 0, 1, . . . . 
The numbers hi and fii are known. Taking r = a in (1.2) and approximating ^f^i by 

under the sign of the integral one gets: 

(1.4) / q{p)fi{p)dp = bu / = 0,1,2,3,... 
Jo 

where 

fi{p) ■= mia,p)ipl"\p)- 

Differentiate (1.2) with respect to r, set r = a, and again replace ipi by tp^^^ under 
the sign of the integral. The result is: 

(1.5) / q{p)si{p)dp = Pi, ^ = 0,1,2,3,... 
Jo 

where 

(1.6) Mp):=lim^ ^l'\p). 

or r=a 

We have replaced ipiip) in (1.2) under the sign of the integral by ipl^\p). Such 
an approximation is to some extent similar to the Born approximation and can be 
justified if q{p) is small or I is large or a is small. Note that this approximation is 
also different from the Born approximation since ipj^^ is defined by formula (2.7) 
and incorporates the information about q{r) on the interval r > a. 

For potentials which are not small and for I not too large such an approximation 
cannot be justified theoretically, but may still lead to acceptable numerical results 
since the error of this approximation is averaged in the process of integration. 

We have derived approximate equations (1.4) and (1.6). From these equations 
one can find an approximation to q{p) numerically. These equations yield a moment 
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problem which has been studied in the literature. In |^ and ||T^, Section 6.2], a 
quasioptimal numerical method is given for solving moment problems with noisy 
data. 

Note that the functions s;(p) differ only by a factor independent of p from the 
functions fi{p) defined in (1.5). This is clear from the definition of these functions 
and formula (A. 12) for gi. Therefore one can use equations (1.4) for the recovery of 
q{r), and equations (1.6) are not used below. From the definition of //(p) it follows 
that the set {fi{p)}o<i<L of these functions is linearly independent for any finite 
positive integer L. Also, one can prove that the moment problem (1.4) has at most 
one solution. Indeed, the corresponding homogeneous problem (1.4), corresponding 
to 6/ = for all Z = 0, 1, 2, ... , has only the trivial solution because the set of 
functions {//(p)}o<;<cx3 is complete in L^(0,a) as follows from the result in |^ (see 
also Q). To see this, note that the function {/;(p)}o<z<oo differs only by a factor 
independent of p from the function uf{p) (the functions ui are defined in (A.l)), 
and the set of these functions is the set of products of solutions to homogeneous 
equation (2.2) (in Section 2) for all I = 0, 1, 2, ... . The set of these products is 
complete in L^(0,a) because the set of functions {uf (p)p~^i1(a)y/(/3)}/=o,i,2,...j 
Q!,/3 e S'^, where 5^ is the unit sphere in R^, is the set of products of solutions 
to the homogeneous Schrodinger equation which is complete in LP'{Ba), as follows 
from the results in |l0| . 

Let us explain the idea of the method discussed in detail in , which is similar 
to the well-known Backus- Gilbert method Fix a natural number L and look 
for an approximation of q{r) of the form 

(1.7) q^[r) ■.^Y.^iui[r) := A^ir, p)q{p)dp 

1=0 -^0 

where the kernel is defined by formula (1.10) below, hi are the known numbers 
given in (1.3), and vi{r) are not known and should be found for any fixed r e [0, a] 
so that 

(1.8) \\qL{r) - q{r)\\ ^ as L ^ oo. 

The norm in (1.9) is L^[0, a] or C[0, a] norm depending on whether q e £^[0, a] or 
q e C[0, a]. Condition (1.9) holds if the sequence of the kernels 

L 

(1-9) Ai(r,p) :=^i.i(r)/i(p) 

1=0 

is a delta-sequence, that is, 

(1.10) AL{r,p)^6{r- p), L ^ +x> 

where d{r — p) is the delta-function. 

For (1.11) to hold, we calculate i^i(r) and pi{r) from the conditions: 



(1.11) / AL{r,p)dp=l 

Jo 

(1.12) r\ALir,p)W-prdp = 



= mm . 



"'0 

The parameter r e [0, a] in (1.12), (1.13) is arbitrary but fixed. Condition (1.12) 
is the normalization condition, condition (1.13) is the optimality condition for the 
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delta-sequence Ai(r, p). It says that ^^(r, p) is concentrated near r — p, that 
is, Al is small outside a small neighborhood of the point p — r. One can take 
7 = 2 in (1.13) for example. The choice of 7 defines the degree of concentration 
of ^^(r, p) near p — r. We have taken \Al\'^ in (1.13) because in this case the 
minimization problem (1.12), (1.13) can be reduced to solving a linear algebraic 
system of equations. 

In order to calculate ^^(r) by formula (1.8) one has to calculate vi{r) and /i;(r) 
for different values of r £ [0, a]. 

The problem (1.12), (1.13) is a problem of minimization of the quadratic form 
(1.13) with respect to the variables i'i{r) (which are considered as numbers for 
a fixed r) under the linear constraint (1.12). Such a problem can be solved, for 



example, by the Lagrange multipliers method |10, Section 6.2]. 

In Section 2 we discuss numerical aspects of the proposed method for solving the 
inverse scattering problem. The idea of the method is also applicable to the inverse 
scattering problem with data given at a fixed I for all A; > 0, or for some values 
of k. We specify in equation (1.2) and give a method for computing the data 
(1.1). In Section 3 a summary of the proposed method is given. In the Appendix 
we have collected all the necessary reference formulas in order to make this paper 
self-contained. 

In numerical results obtained by the proposed approximate inversion method 
are given. 

2. Numerical aspects 

2.1. Calculating -0oi and the data (1.1). We start with a method for calculating 
the data (1.1) from the original data {5i, / = 0, 1, 2, . . . } and {q{r), r > a}. 
The following integral equation is convenient for finding the data (1.1): 

/"OO 

(2.1) ipi{r) ^ -iiJoiir) - ^i{r, p)q{p)ipi{p)dp, r>a 

J r 

where the Green function ^i{r, p) is defined in (A. 10) and ipQi is defined by formula 
(2.6) below. Equation (2.1) is a Volterra integral equation. It can be solved stably 
and numerically efficiently by iterations on the interval r > a where q(r) is known. 
The solution to equation (2.1) is used in equation (2.7) to determine V'^'^'- 
The function tpoi is the unique solution to the equation 

(2.2) Voi + fcVoi - ^^^^V^oi =0, r > 
which has the same asymptotics as "0; as r — > -l-oo: 

(2.3) V/ - e'^' sin (fcr - y + 5/) , r ^ 00. 

This formula is derived in Appendix (formula (A. 29)). Let us derive a formula for 
ipoi. If ipai has the asymptotics (2.3), then 

(2.4) xpoi ^ ciui{kr) + C2Vi{kr) ci sin {kr -—)- C2 cos [kr r ^ 00 

where ui and vi are defined in (A.l) and ci,C2 are some constants. 
From (2.4) and the asymptotics (2.3) for tpoi, it follows that 

(2.5) ci = e'^' cos(,5j), C2 = -e"' sin((5j). 



INVERSE SCATTERING PROBLEM 5 

Thus we obtain an explicit formula for ipoi- 

(2.6) ^oiikr) = e"'' cos{Si) ui{kr) - e'^' sm{Si) vi{kr). 

2.2. Calculating V'!''^''- Let us find ipl^'' in equation (1.2). We start with the in- 
tegral equation (A. 14) for the physical wave function ipi and rewrite (A. 14) as 
equation (1.2) with 

(2.7) =i,l"\r,k) ■.= ui(kr) - H gi{r, p)q{p)Mp, k)dp. 

J a 

Formula (2.7) solves the problem of finding -^j-"-' in equation (1.2). Indeed, in 
Section 2.1 we have shown how to calculate ipi{p,k) for p > a, so that the right- 
hand side of (2.7) is now a known function of r. This completes the description of 
the numerical aspects of the proposed inversion method. 



3. Summary of the proposed method 
Let us summarize the method: 

1. Given the data {Si, Z = 0, 1, 2, . . . ; q{r), r > a} one calculates "00/ by formula 
(2.6), solves equation (2.1) by iterations for r > a, and obtains the data (1.1). 

2. Given the data (1.1), one calculates bi by formulas (1.3), then chooses an in- 
teger L and solves the moment problem (1.4) for < I < L. The approximate 
solution qhi'i') is given by formula (1.8) in which vi{r) are found by solving 
the optimization problem (1.12), (1.13). 

In conclusion we make several remarks. 

Remark 1. This method can be formulated also for the case when the data hi of 
the moment problem are known with some random errors (see jl^, - Section 6.2] 
for details). This idea can also be tried for a numerical inversion of fixed-energy 
scattering data. 

Remark 2. This inverse scattering problem with fixed-energy data is highly ill- 
posed if k is not very large. Specific estimates of the conditional stability of the 
solution to this problem were obtained in pl| and illustrated in [Q (sec also pO[|). 
Although it was proved in |^ (see also 0-[^) that the 3D inverse scattering problem 
with fixed-energy scattering data has at most one solution in the class of compactly 
supported potentials, but the estimates of the stability of this solution with respect 
to perturbation of the scattering data obtained in |1^ indicate that the stability is 
very weak: theoretically one has to have the scattering data known with very high 
accuracy in order to recover the potential with a moderate accuracy. Of course the 
estimates in ||ll|] cover the worst possible case, and in practice the recovery may be 
much better than the error estimates from |]ll| guarantee. 

In a recent paper Q an example is given of two quite different piecewise-constant, 
real-valued, compactly supported potentials which have practically the same phase 
shifts at a fixed energy for all I (the phase shifts differ by a quantity of order 10~^). 

This result shows that it is necessary to have some a priori information about 
the unknown potential q{r) in order to be able to recover it with some accuracy. 

Remark 3. In the Born approximation, inversion of fixed-energy scattering data 
is reduced to solving an integral equation of the form: 

(3.1) / q{x)e-'^-^dx = f{i), |e|<2fc 
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where fc > is a fixed given number, x £ R^, and a > is the radius of a ball Ba 
outside of which q{x) — 0. If one knows /(^) for all ^ G R^, one can find q{x) by 
taking the inverse Fourier transform of /(^). If /(^) is known only for |^| < 2fc, 
one can still recover a compactly supported qix) uniquely and analytically with an 
arbitrary accuracy, if /(^) is known exactly for |^| < 2k (see Section 6.1], where 
analytic inversion formulas for (3.1) are derived). 

Remark 4- A numerical implementation of the parameter- fitting procedure of 
R. G. Newton requires a priori information about the potential q, in particular, the 
knowledge of g for r > a, see 

4. Appendix 

Equation (2.2) has linearly independent solutions uiikr) and vi{kr), called Ric- 
cati-Bessel functions: 

/ 7Tr / 7Tr 

(A.l) ui{r) J Y J;+i/2(r) := rji{r), vi{r) = J y^i+iM^) rni{r) 

where Jiy{r) and N^{r) are the Bessel and Neumann functions regular and irregular 
at the origin respectively, ji and ni are the spherical Bessel and Neumann functions. 
One has 

In In 
(A. 2) vi{r) ^ - cos{r - Mi(r) ~ sin(r - y), r — » cc 

f2; - IV r'+^ 
(A.3) „,(,)^-^_J::, ^,(,) ^ , ^ o. 

The regular solution to (2.2), denoted by ipQi{kr) and defined by the condition 

{2l + l)\V 



(A.4) mikr)^———, r^O 



is given by the formula 

(A.5) ^o^ikr) = ^ 

which follows from (A.3) and (A.4). 

The Jost solution foi to (2.2) defined by the condition 

(A.6) /oi(fcr-) -e^"^ r ^ cx) 

is given by the formula 

(A.7) Mkr) - e'('+i)^/2^,(fcr) - e''^/\iikr) = ie''^/^{ui + ivi) 



as follows from (A. 2). 
The Wronskian is 



(A.8) Foiik) := W[foi,m] = foi{kr) — ipoi{kr) - m{kr) — foi{kr) = 

as follows from (A.3), (A.4), (A.7), and the fact that this Wronskian does not 
depend on r, so it can be calculated at r = 0. 

In equation (2.1) we use the Green function ^i{r, p) satisfying the equation 

A.9) — — + k^i{r,p) 2 — iiif, P) ^ -S{r - p) 

ar'' r'^ 
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(A. 10) ii{r,p) 



where 5{r — p) is the delta- function, and vanishing for p < r: 

_ / e-'^^/''k^[-ipoi{kp)foi{kr) + foi{kp)<poi{kr)], p > r 

0, p<r. 

Indeed, solves (A. 9) for r 7^ p since, as a function of r, it is a linear combination 
of solutions to the homogeneous equation (A. 9). Moreover, 

(A.ll) 



Voi{kp) J- foi{kp) 



dr dr 



r=p 



Umi^^^ =e-^'-/'k'lun 

dr r=p-o 

by (A.8). Thus (A.9) follows. 
The Green function 

(A 12) oAr n\ - / ^o7'W</'oj(M/o/(M> r>p 

^ ' ^'^''''^-\ F,-,\k)<poi{kr)foi{kp), r<p 

(where FQi{k) is defined in (A.8), Fqi = e'''^/^fc~') solves (A.9) and satisfies the 
radiation condition at infinity: 

(A. 13) — ikgi as r — »■ +00. 

or 

Indeed, (A. 13) follows from (A. 12) and (A. 6), and to verify that gi solves equation 
(A.9) one argues in the same way as was done after formula (A. 10). 
The integral equation for the physical wave function ipi (r, k) is: 



(A.14) i)i{r, k) = ui{kr) - gi{r, p)q{p)i)i{p, k)dp. 

Jo 

This equation follows from the standard three-dimensional equation: 
(A.15) 

il){x,a,k) =il)o{x,a,k) - I — -q{p)tp{y,a,k)dy, p=|y|, r=|x| 

Jr3 47r|a; - y\ 

where 'ipo{x, ct, k) := e'^"'^, and a G 5^ is a given unit vector. 
If one writes 

(A.16) e^ka:^y^^iMkr)_ o^y^ x° := 

^ k r \x\ 

(A.17) ^ = Ef^'^^K^°)W 

(A.18) = V ^FK-") W) 

where Yi := Yjm, —l<m<l, are the orthonormal spherical harmonics, and inserts 
(A.16)-(A.18) into (A.15), then one gets (A.14). In (A.16)-(A.18) the summation 
with respect to Z, Z = 0, 1, 2, 3, . . . , includes also the summation with respect to m, 
-l<m<l. 

Let r ^ +00 in (A.14). Then (A.2), (A.5), (A.6), (A.8), and (A.12) imply 

„j7r(i+l)/2 

(A.19) [e-''='- - e'^'e'^'-Si], r ^ +00 
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where 

(A.20) Si = l + — / ui{kp)q{p)iji{p, k)dp. 

Jo 

The scattering ampUtude A{a',a,k) defined by the formula 



oikr 



(A. 21) ij; = escpUka ■ x) + A{a' ,a,k) l-o(-), r := |a;| — »■ oo, — = a' 

r r \x\ 

can be written as 

(A.22) A{a', a, fc) = --1 / e-*'=«'-^(z(p)V(y, a, k)dy. 

If gr = q{p), p := \y\, then 

OO 

(A.23) A{a',a,k) = Y,Mk)Yi{a')Wa) 

1=0 

where 

47r /"^ 

(A.24) Ai{k) = —^j^ uiikp)q{p)^i{p,k)dp 

as one gets after substituting the complex conjugate of (A. 16), (A. 17), and (A.23) 
into (A.22). From (A.20) and (A.24) one gets 

(A.25) Si = l- ^Ai 

which is the fundamental relation between the iS-matrix and the scattering ampli- 
tude: 

(A.26) S = I- —A. 

Here S, /, and A are operators in i^(S'^), / is the identity operator. Since the 
operator S is unitary in LP'{S'^) and Si are the eigenvalues of S in the eigcnbasis of 
spherical harmonics in the case of spherically symmetric potentials, one has \Si\ = 1, 
so, for some real number one writes 

(A.27) Si = e^'^K 

The numbers i5; normalized by the condition 

(A.28) (5i(fc)^0 as A; ^ oo 

arc called the phase shifts. They are in one-to-one correspondence with the numbers 
Si if (A.28) is assumed. Wo do assume (A.28). 
Formula (A. 19) can be rewritten as 

(A.29) 

^i(kr-l-K 11-Vii) _ ^-i(kr-l-K/2+5i) 

ipi ~ e"'' — = e*''' sin (kr - — + 6i) , r ^ oo. 

Formula (A.29) was used in Section 2 (see formula (2.3)). 
Note that (A.25) and (A.27) imply 

(A.30) Ai = ^e'''sm{Si). 

This formula and (A.23) show that if q{x) = q{\x\), then the information which 
is given by the knowledge of the fixed-energy phase shifts {di, I = 0, 1,2, . . .} is 
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equivalent to the information which is given by the knowledge of the fixed-energy 
scattering amplitude A{a',a) where a' and a run through the whole of the unit 
sphere 5^. By Ramm's uniqueness theorem it follows that if the fixed-energy 
phase shifts corresponding to a compactly supported potential q{r) are known for 
I — 0,1,2,3, ... , then the potential q{r) is uniquely defined, and therefore the phase 
shifts di are uniquely defined for all I on the complex plane by their values for I 
running through the integers only: 1^0,1,2,.... 
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